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Galactic halo neutralinos (%) captured by the Sun or Earth produce high-energy 
neutrinos as end-products of various annihilation modes. These neutrinos can travel 
from the Sun or Earth cores to the neighborhood of underground detectors ("neu- 
trino telescopes" ) where they can interact and produce upwardly-moving muons. We 
compute these muon fluxes in the context of the minimal 577(5) supergravity model, 
and the no-scale and dilaton 577(5) x £7(1) supergravity models. At present, with the 
Kamiokande 90% C.L. upper limits on the flux, only a small fraction of the parameter 
space of the SU(5) x £7(1) models is accessible for m x ~ m-F e , which in turn implies 
constraints for the lightest chargino mass around 100 GeV for a range of tan f3 values. 
We also delineate the regions of parameter space that would be accessible with the 
improvements of experimental sensitivity expected in the near future at Gran Sasso, 
Super-Kamiokande, and other facilities such as DUMAND and AMANDA, currently 
under construction. We conclude that if neutralinos are present in the halo, then this 
technique can be used to eventually explore more than half of the allowed parameter 
space of these specific models, and more generally of a large class of supergravity 
models, in many ways surpassing the reach of traditional collider experiments. 
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1 Introduction 



Even though there appear to be several indications that low-energy supersymme- 
try is indeed a symmetry of Nature JI[, no more converts should be expected until 
an actual sparticle is observed experimentally or unequivocal indirect experimental 
evidence becomes compelling. For either discovery mode to be efficient, one must 
provide as accurate experimental predictions as possible. In supersymmetric models 
where the number of parameters is large, as in the minimal supersymmetric standard 
model (MSSM), easily falsifiable experimental predictions are hard to obtain since 
the various parameters usually allow a myriad of possibilities. One basic step forward 
is to study specific supersymmetric models wherein the assumptions made are well 
motivated by physics at very high energies, such as grand unification, supergravity, 
and superstrings. In this way one can hope to test general theoretical frameworks 
experimentally, as opposed to just ruling out small regions of a many-dimensional 
parameter space. This line of thought has been in existence for almost as long as 
supersymmetric phenomenology has. However, it has only been since the firm estab- 
lishment of the standard model at LEP that real candidates for the theory beyond 
the standard model have emerged from a morass of pre-LEP challengers. 

The aim of this paper is to study a very promising indirect experimental proce- 
dure to discover supersymmetry, in the context of a class of SU (5) and SU (5) x U (1) 
supergravity models, via interactions of upwardly-moving muons in underground de- 
tectors. As we will show, this method of experimental exploration is quite competi- 
tive with other indirect probes which have been recently investigated in the context 
of this same class of models, namely supersymmetric contributions to the one-loop 
electroweak LEP observables ||, to the rare radiative b — > 57 decay ||, and to the 
anomalous magnetic moment of the muon [|J. 

The basic idea is that the neutralinos \ (lightest linear combination of the su- 
perpartners of the photon, Z-boson, and neutral Higgs bosons), which are weakly 
interacting massive particles (WIMPs), are assumed to make up the dark matter in 
the galactic halo, and can be gravitationally captured by the Sun or Earth || 
after losing a substantial amount of energy through elastic collisions with nuclei. 
The neutralinos captured in the Sun or Earth cores annihilate into all possible ordi- 
nary particles, and the cascade decays of these particles as well as their interactions 
with the solar or terrestrial media produce high-energy neutrinos as one of several 
end-products. These neutrinos can then travel from the Sun or Earth cores to the 
vicinity of underground detectors, and interact with the rock underneath producing 
detectable upwardly-moving muons Such detectors are rightfully called "neutrino 
telescopes" , and the possibility of indirectly detecting various WIMP candidates has 
been considered in the past by several authors [0] . More recent analyses can be found 
inRefs. §^0,0,0,0,0. 



^dn this paper, we do not consider the less promising "contained events" in which neutrino 
interactions take place within the detector, since the event rate for such events is proportional to 
E u , as opposed to E„ for the upwardly- moving muon events. 
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Making use of the techniques in the literature [|| [10], IT]], our current work focuses 
on the predictions for the upwardly-moving muon event rates (fluxes) in two dis- 
tinct supersymmetric models, namely, the minimal 577(5) supergravity model and the 
string-inspired S77(5) x U(l) supergravity model. In these models, because of the un- 
derlying structure of supergravity, the low-energy couplings and masses are completely 
determined via renormalization group equations (RGEs) in terms of only five param- 
eters: the three universal soft-supersymmetry-breaking parameters (myz, tuq, A), the 
top-quark mass (m t ), and the ratio of Higgs vacuum expectation values (tan/5). The 
restriction to just five parameters is possible because of the radiative electroweak 
symmetry breaking mechanism, which allows to determine some further parameters, 
like the magnitude of the Higgs mixing term [i (but not its sign). This mechanism 
imposes further constraints on the model parameters (e.g., it requires tan/3 > 1) and 



involves the value of m t in a fundamental way p5| . These features make it possible 
to incorporate all currently available experimental constraints, both direct and indi- 
rect, as well as some theoretical consistency conditions into a complete determination 
of the allowed parameter space of these models. We calculate the upwardly- moving 
muon fluxes induced by the neutrinos from the Sun and Earth in the still-allowed pa- 
rameter space of these models, and compare them with the currently most stringent 
90% CL. experimental upper bounds, obtained at Kamiokande, for neutrinos from 
the Sun [IIJ and Earth Jl!| respectively, i.e. 

Tsun < 6.6 x lO^cm-V 1 = 2.08 x l(T 2 iir 2 yr _ \ (1) 
r E arth < 4.0 x lO-^cm-V 1 = 1.26 x lO^m-V -1 - (2) 

Aiming at the next generation of underground experimental facilities, such as MACRO 



and other detectors at the Gran Sasso Laboratory [17], Super-Kamiokande [18], DU- 



MAND, and AMANDA [IP], where improvements in sensitivity by a factor of 2-100 
are expected, we also delineate the region of the parameter space of these models that 
would become accessible with an improvement of experimental sensitivity by modest 
factors of two and twelve. 

This paper is organized as follows. In Sec. § we describe the minimal SU(5) 
and SU(5) x £7(1) supergravity models. In Sec. ^| we outline the calculation of the 
capture rates at the Sun and Earth, while in Sec. |4] we outline the calculation of the 
corresponding detection rates. In Sec. |5| we present the results of our computations 
for the two supergravity models described in Sec. |2|. Finally, we conclude in Sec. |] 
with some comments. 



2 The Supergravity Models 

We work in the context of two supergravity models based on the gauge groups SU(5) 



23 and SU(5) x 17(1) ("flipped S77(5)") 0. In these two models, spontaneous 
supersymmetry breaking takes place in the "hidden" sector, which manifests itself 
in the "observable" sector Lagrangian as a set of soft-supersymmetry-breaking terms 
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that are universal at the respective unification scales of the models. The low-energy 
particle content of these two models is the same as that of the MSSM. However, 
mainly because of the different gauge group structure, these two models have rather 
different low-energy phenomenologies, which we now briefly describe in turn. 

In the minimal SU(5) supergravity model, the unification of the standard model 



gauge couplings takes place at a scale My ~ 10 GeV p2[, as a combined result of the 



SU (5) gauge group and the minimal matter content. This property of the model has 
been shown to hold even after the most general low-energy and high-energy thresh- 
old corrections have been incorporated, as well as the two-loop expressions for the 
RGEs [21|. The unified symmetry implies the existence of additional heavy particles 
with masses of order M\j. In particular, the Higgs sector includes at least the 24 
representation to break SU(5) down to SU(3) x 577(2) x £7(1), as well as the 5, 5 
superfields, whose scalar doublet components (-££2) are responsible for the electroweak 
gauge symmetry breaking and should be kept light, while the scalar triplet compo- 
nents (-££3) must be heavy to avoid fast proton decay through dimension-six operators. 
It also follows that the Yukawa couplings of the bottom-quark and tau-lepton must 
be unified at Mjy, which when evolved to low energies entails a constraint on the rela- 
tion between the ratio of Higgs vacuum expectation values (tan f3) and the top-quark 
mass {nit). Perhaps the most distinguishing feature in the minimal SU(5) model has 
to do with the genuinely supersymmetric proton decay via dimension-five operators, 
mediated by the exchange of superpartners of the heavy Higgs triplets (-££3). Even 
for sufficiently large values of M^ 3 , it is necessary to tune the sparticle spectrum 
so that this type of proton decay remains at an acceptable level |23|| . This usually 
requires light charginos and neutralinos, while squarks and sleptons should be heavy. 
We study a generic supersymmetry breaking scenario characterized by the universal 
parameters (mi/ 2 ,m , A). It turns out that, because of the correlation between the 
relic density of the lightest neutralino and the soft supersymmetry breaking patterns 
|[24|| , the cosmological constraint together with the proton decay constraint dramati- 
cally reduce the parameter space of the minimal SU(5) model, and what remains are 
essentially points corresponding to the resonances in the neutralino pair annihilation 
cross section where m x ~ \Mz, \vrih f25fl . In addition, the ratio £0 = mo/ m i/2 must 



be significantly larger than unity, and tan/5 must not be too large (tan/3 ^3 — 5) 



[p3|, p6| . The ensuing simple relations among the different sparticle masses in this 
model can be summarized as: 

Higgs : 60 GeV < m h < 125 GeV, 
X1.2, Xi ■ 2m ° ~ m x % ~ «\± ~ 0.3m 5 , 



Minimal SU(5) 



Squarks : rriq ~ \m g \J§ + £0, 

Sleptons : m e - fi « \rrigJ0.15 + £q, ps \m g Jo.5 + £q 



(3) 

where £ ^ 3-6 depending on the value of the triplet higgsino mass used (M^ a < (3- 
10)Mjj). The actual five- dimensional parameter space to be explored in this paper 
for the minimal SU(5) model has been obtained in Ref. p6f . 

In the string-inspired SU(5) x £7(1) supergravity model, the preferred unification 
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scale is Mjj ~ 10 18 GeV, as expected in the context of string theory P?fl , which has 
been realized by the effects of extra vector-like matter fields at some intermediate 
scales |[28|| . Contrary to the minimal SU(5) model, the Higgs bosons needed to break 
the SU(5) x U(l) symmetry belong to the 10 and 10 representations, which can be 
easily accommodated by the simplest string models (with Kac-Moody level k = 1), 
while the adjoint representation (24) may only appear in more complicated construc- 
tions. Because of the SU (5) x U{1) gauge group structure, the doublet-triplet splitting 
of the pentaplet Higgs superfields occurs in a very simple way, and the dimension-five 
proton decay operators are automatically suppressed. In the case of the SU (5) x U (1) 
supergravity model, we study two string-inspired supersymmetry breaking scenarios: 
(i) the no-scale model J29|, where m = A = f3(J, and (ii) the dilaton model f31~H 



where m = A = — mi/ 2 fl32] . Therefore, the SU(5) x U{1) model (either sce- 

nario) is quite predictive since it depends on only three parameters: m t , tan/3, mi/ 2 . 



Also, in these two scenarios the cosmological constraint is automatically satisfied |24] 
The mass relations in this model are: 



SU(5) x U(l) : < 



Higgs : 60 GeV < m h < 125 GeV, 

Xi,2, xf ■ 2m x o ss m x o « m % ± « 0.28m 5 , 

Squarks : rrig m rrig, 

Sleptons : m~ eR » 0.18(0.33)m 5 , mg L « 0.30(0.41)m g , 



for the no-scale (dilaton) case. Note that the squark and slepton mass relations 
in Eq. (^) do not reduce to those in Eq. (TJ), for values of £ — (no-scale) and 

= l/y/3 (dilaton), because of slight changes in the running of the scalar masses 
down to low energies from the different starting value of My. The allowed three- 
dimensional parameter space of this model has been determined in Refs. [^, for 
the no-scale and dilaton cases respectively. 

We now comment on some other features about these two models which are par- 
ticularly relevant to our current work. First, for m g ^ 1 TeV the lightest neutralino 



can not be a pure higgsino |33[]. In fact, we have found that for almost all the points 
in the allowed parameter spaces of both models, the lightest neutralino is a "mixed" 
state with substantial gaugino and higgsino components. In the MSSM, such "mixed" 



neutralinos normally cannot even account for the halo dark matter p4|, because the 



efficient annihilation via Higgs boson exchange renders the relic density rather small 
(fi x /ig < 0.05). This situation is improved in the two supergravity models, as a result 
of the one-loop radiative corrections to the masses of Higgs bosons which we have 
included in our analysis. Since the one-loop corrected Higgs masses are normally 
larger than the tree-level values, the overall effect is an enhancement of the relic den- 
sity for "mixed" neutralinos. Therefore, the lightest neutralinos in these two models, 
although mostly "mixed" states, are still good candidates for the major component 
of the galactic halo. Such neutralinos are mainly captured by the Sun and Earth 
through their coherent (spin-independent) scattering off nuclei due to the exchange 
of Higgs bosons. Furthermore, in these two models the neutralinos typically have 
masses in the 20-150 GeV range, so their capture by the Earth is expected to be 
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enhanced when the neutralino mass closely matches that of the abundant elements 
in the Earth's core (Fe) and mantle (Si and Mg), while the same effect is irrelevant 
for the capture by the Sun 0. We will discuss the implications of these features in 
Sec. § 



3 The Capture Rate 

In order to calculate the expected rate of neutrino production due to neutralino 
annihilation, it is necessary to first evaluate the rates at which the neutralinos are 
captured in the Sun and Earth. Following the early work of Press and Spergel ||, 
the capture of WIMPs by a massive body was studied extensively by Gould . In 
this paper, we make use of Gould's formula, and follow a procedure similar to that 
of Refs. [11]] in calculating the capture rate. 

From Eq. (A10) of the second paper in Ref. ||, the capture rate of a neutralino 
of mass m x by the Sun or Earth can be written as 

G= (2y^PA E A <7 , Xj , (5 ) 

where Mb is the mass of the Sun or Earth, p x and v x are the local neutralino density 
and rms velocity in the halo respectively, cij is the elastic scattering cross-section of 
the neutralino with the nucleus of element i with mass rrii, fi is the mass fraction of 
element i, and A; is a kinematic factor which accounts for several important effects: 
(1) the motion of the Sun or Earth relative to Galactic center; (2) the suppression 
due to the mismatching of m x and mf, (3) the loss of coherence in the interaction due 
to the finite size of the nucleus (see Ref. || for details). 

In the summation in Eq. (|5]), we only include the ten most abundant elements 
for the Sun or Earth respectively, and use the mass fraction fi of these elements 
as listed in Table A.l of Ref. [TO]. We choose v x = 300km value within the 



allowed range of the characteristic velocity of halo dark matter particles. To take into 
account the effect of the actual neutralino relic density, we follow the conservative 
approach of Ref. |TOJ for the local neutralino density p x : (a) p x = ph = 0.3 GeV/cm 3 , 
if Q x h 2 > 0.05; while (b) p x = {Q x h 2 /0.05)p h , if Q x h 2 <, 0.05. As for a u the 
dominant contribution is the coherent interaction due to the exchange of two CP- 
even Higgs bosons h and H and squarks, and we use the expressions (A10) and 
(All) of Ref. |lTji to compute the spin-independent cross section for all the elements 



included. In addition, for capture by the Sun, we also evaluate the spin-dependent 
cross section due to both Z-boson exchange and squark exchange for the scattering 



from hydrogen according to Eq. (A5) (EMC model case) of Ref. [0]]. It should be 
noted that in all these expressions the squarks were assumed to be degenerate. In the 
two supergravity models that we consider here this need not be the case, although for 



2 We have corrected a sign error for the H\X coupling in (A10) of Ref. 1 1 1 
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most of the parameter space this is a fairly good approximation. Hence, we simply 
use the average squark mass rrig in this part of the calculation. 

The kinematic factor X, in Eq. (|5|), can be most accurately evaluated once the 
detailed knowledge of the mass density profile as well as the local escape velocity 
profile are specified for all the elements. In practice, this can be done by performing 
a numerical integration with the physical inputs provided by standard solar model 
or some sort of Earth model. Instead of performing such an involved calculation, we 
approximate the integral for each element by the value of the integrand obtained with 
the average effective gravitational "potential energy" fa times the integral volume. 



The values of fa are taken from Table A.l of Ref. 10 



4 The Detection Rate 

We next describe the procedure employed by us to calculate the detection rate of 
upwardly-moving muons, resulting from the particle production and interaction sub- 
sequent to the capture and annihilation of neutralinos in the two supergravity models. 
The annihilation process normally reaches equilibrium with the capture process on a 
time scale much shorter than the age of the Sun or Earth. We assume this is the case, 
so that the neutralino annihilation rate equals half of the capture rate. The detection 
rate for neutrino-induced upwardly-moving muon events is then given by 

" C EA^/W EldE v . (6) 



8nR 2 iF , \—v, iF 

In Eq. Di is a constant, R is the distance between the detector and the Sun or the 
center of the Earth, and [dN / dE u )ip is the differential energy spectrum of neutrino 
type i as it emerges at the surface of the Sun or Earth due to the annihilation of 
neutralinos in the core of the Sun or Earth into final state F with a branching ratio 
Bp. It should be noted that in Eq. (|6|) that i is summed over muon neutrinos and 
anti-neutrinos, and that F is summed over final states that contribute to the high- 
energy neutrinos. The only relevant fermion pair final states are rf , cc, bb, and (for 
the SU(5) x U(l) model) tt when m x > m t . The lighter fermions do not produce 
high-energy neutrinos since they are stopped by the solar or terrestrial media before 
they can decay ||. 

The branching ratio Bp can be easily calculated as the relative magnitude of the 
thermal-averaged product of annihilation cross section into final state F (o>) with the 
M0ller velocity vm- Since the core temperatures of the Sun and Earth are very low 
compared with the neutralino mass (T Sun ~ 1.34 x 1CT 6 GeV; T Earth ~ 4.31 x 1CT 10 
GeV), only the s-wave contributions are relevant, hence it is enough here to use the 
usual thermal average expansion up to zero-order of T/m x {vm — ► limit), i.e. 

D (ctfVm) a F 

Bp = = . (7) 
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In Eq. (0), all kinematically allowed final states contribute to Otot- Besides all the 
fermion pair final states, we have also included boson pair final states WW, ZZ and 
hA in our calculation of Bp. Due to the parameter space constraints, these channels 
are not open for the minimal SU(5i) model. But WW and ZZ channels are generally 
open in the SU(5) x U{1) model, and the hA channel also opens up for large values 
of tan/3 in the dilaton case. We should also remark that the annihilation channel 
into lightest CP-even higgs pair hh is always allowed kinematically in some portion 
of the parameter space for both supergravity models we consider, but since its s-wave 
contribution vanishes, we do not include it in the calculation of Bp. However, this 
channel is taken into account in the calculation of the neutralino relic density, which 
does affect the capture rate through the scaling of local density p x when Vt x h1 < 0.05 
(see Sec. [3|). In addition, we have kept all the nonvanishing interference terms in the 
evaluation of a tot and dp. 

The calculation of the neutrino differential energy spectrum is somewhat involved, 
since it requires a reasonably accurate tracking of the cascade of the particles which 
result from neutralino annihilation into each of the final state F. This involves the 
decay and hadronization of the various annihilation products and their interactions 
with the media of the Sun or the Earth's cores. In addition, at high energies, neutri- 
nos interact with and are absorbed by solar matter, a fact that affects the spectrum. 
In Ref. ||, Ritz and Seckel rendered this calculation tractable by their adaptation 
of the Lund Monte Carlo for this purpose. Subsequently, analytic approximations 
to the Monte Carlo procedure outlined in their paper were refined and employed by 
Kamionkowski to calculate the neutrino energy spectra from neutralino annihilation 
for the MSSM in Ref. [IT]. The procedure involved is described below for complete- 



ness. 

Since the probability for producing an underground muon that traverses the de- 
tector is proportional to the square of the neutrino energy, the primary quantity of 
interest is the second moment (Nz 2 )m x , defined as 

hp — 

m 




Once the second moments are obtained, the detection rate for the neutrino-induced 
upwardly-moving muon events (0) may be conveniently written as 

r = ^(^ T )(^) 2 E^E^(^ 2 ) iF m-V 1 , (9) 

where Kb = 1.27 x 10 -29 (7.11 x 10~ 21 ) for neutrinos from the Sun (Earth), the 
scattering coefficients <2j = 6.8 (3.1) for neutrinos (anti-neutrinos), while the muon 
range constants bi = 0.51 (0.67) for neutrinos (anti- neutrinos). 



The approximations of Refs. ||, 11] consist in obtaining expressions for the second 
moment without detailed knowledge of the functional form of the differential energy 
spectra. We now list these approximate expressions for (Nz 2 ) iF as follows: 
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(i) Fermion pair final states (rT,cc,bb,tt). The simplest case is that of fermions 
injected into the core of the Earth, since all interactions are negligible. We use 



(Nz 2 ) = -(N)(y 2 ) 



f 



(10) 



Here (N) and (y 2 ) are the rest frame yield and second moments, while (zf) is the 
second moment of the fragmentation function, obtained from Table 2 and Table 3 
respectively in ||. rrif and E in are the mass and energy of the injected fermion. 

Fermions injected into the core of the Sun undergo interactions and decay before 
final state neutrinos emerge at the surface. These are analytically approximated 



reliably for fermion injection energies relevant to our situation by ['I I 



(Nz 2 ) 



-\ 2 



axoe x ° I dx 

XQ X 



where x = 155/ E in (275/ E in ) and a = 0.056 (0.052) for neutrinos (anti- neutrinos) 
from c quarks, while x = 185/ E in {275/ E in ) and a = 0.086 (0.082) for neutrinos (anti- 
neutrinos) from b quarks. For the top quark we have considered three different masses: 
130, 150 and 170 GeV in the SU(5) x £7(1) model. In this case the approximations 



are less reliable [pTTj , but better at energies below the 100 GeV scale than above 
it. The appropriate values here are a = 0.18 (0.14) and x = 110/E in (380/E in ) for 
neutrinos (anti-neutrinos). In the above (and in what follows), E in and other energies 
are taken in GeV when obtaining numerical values. The above expressions include 
the hadronization of b and c and t quarks in the solar medium. The r lepton decays 
almost instantly, and does not hadronize, so must be treated differently The second 
moment for the r is approximated by 



(Nz 2 ) 



ae 



-Ei n /E a 



(12) 



where a = 0.0204 (0.0223) and E a = 476 (599) for neutrino (anti-neutrino) produc- 
tion. 

(ii) Gauge boson pair final states {WW, ZZ). The production of and by 
W, Z boson pair and Higgs boson pair is related to that by fermions, since the bosons 
primarily decay into fermions, hence the subsequent cascade and hadronization re- 
mains the same. For the annihilation final state WW, the W boson with velocity f3 
and energy E will decay into fiu^ with a fractional width T(W — > fiP/i). The second 



moment for this annihilation final state for the neutrinos from the Sun is [11 



(Nz 2 ) 



T(W -> /xz/ M ) 2 + 2E inn {l + on) + E 2 n r 2 ai {l + a 



(3E 3 



airf{a 2 - 1)(1 + E in r i ) 1 +t 



-, E in =E{X-P)/2 



J E m =E{l+0)/2 



(13) 



Here n = 1.01 x 10" 3 GeV" 1 (3.8 x 10" 4 GeV 



r-l\ 



a, 



5.1 (9.0) for neutrinos (anti- 



neutrinos). For the neutrinos from the Earth, the second moments is approximated 
by a simpler expression 



(Nz 2 ) =T{W ^){3 + [3 2 )/12 



(14) 
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due to the absence of interactions. The expressions for the ZZ final state can be 



similarly obtained [11 



(iii) Higgs boson pair final state (hA). For the annihilation final state hA, each 
Higgs boson S (= h, A) with velocity (5s and energy Es will decay into fermion pair 
// with a branching ratio BR(S — > //). Assuming such decays are isotropic, the 
second moment of this final state is 

2BR { s -. fj) i^mn 



Here Ej is the fermion energy, and (Nz 2 )j is the second moment for the fermion pair 
// (rf,cc,bb,ti) as given in Eqs. (PH)-fll^). 



5 Results and Discussion 

For each point in the parameter spaces of the two models described in Sec. [2|, namely 
the minimal SU(5) and the SU(5) x U(l) supergravity models (both the no-scale 
and dilaton cases), we have determined the relic abundance of neutralinos and then 
computed the capture rate in the Sun and Earth (as described in Sec. |3|) and the 
resulting upwardly-moving muon detection rate (as described in Sec. [|). In Figs. |l] 
and 0, the predicted capture and detection rates in the minimal S77(5) supergravity 
model are shown, based on the assumption that the mass of the triplet higgsino, 
which mediates dimension-five proton decay, obeys M# 3 < 3M^. We have redone the 
calculation relaxing this assumption to Mg 3 < lOMu, in which case the results for 
the muon fluxes remain qualitatively the same, except that the parameter space is 
opened up somewhat. The dashed lines in Fig. |2] represent the current Kamiokande 
90% C.L. upper limits (1) and (2). Similarly, the predictions of the SU(5) x U(l) 
supergravity model are presented in Figs. |3| and |] for the no-scale scenario, and in 
Figs. |5] and ^| for the dilaton scenario, again along with the Kamiokande upper limits 
(dashed lines). In Figs. we have taken the representative value of m t = 150 GeV. 
Similar results are obtained for other values of mt- 

Several comments on these figures are in order. First, the kinematic enhancement 
of the capture rate by the Earth manifests itself in all figures as the big peaks near the 
Fe mass (mp e = 52.0 GeV), as well as the smaller peaks around Si mass (msi = 26.2 
GeV). Second, there is a severe depletion of the rates near m x = \Mz in Figs. [3|-§, 
which is due to the decrease in the neutralino relic density. In the case of Earth 
capture, this effect is largely compensated by the enhancement near the Fe mass. As 
mentioned in Sec. in our procedure, the relic density affects the local neutralino 
density p x only if £\/ig < 0.05, while in the minimal SU{5) model this almost never 
happens, therefore, the effect of the Z-pole is not very evident in Figs. [1] and |[ Also, 
in Figs. the various dotted curves correspond to different values of tan (3, starting 
from the bottom curve with tan (3 = 2, and increasing in steps of two. These curves 
clearly show that the capture and detection rates increase with increasing tan /3, since 



9 



the dominant piece of the coherent neutralino-nucleon scattering cross section via the 
exchange of the lightest Higgs boson h is proportional to (l+tan 2 /3). The capture rate 
decreases with increasing m x , since the scattering cross section falls off as m^ 4 and 
rrih increases with increasing m x . It is expected that the detection rate in general also 
decreases for large value of m x , since it is proportional to the capture rate. However, 
the opening of new annihilation channels, such as the WW, ZZ and HA channels in 
the SU(5) x U(l) model, could have two compensating effects on the detection rate: 
(a) the presence of a new channel to produce high-energy neutrinos which leads to an 
enhancement of the detection rate; and (b) the decrease of the branching ratios for the 
fermion pair channels, which makes the neutrino yield from tt, cc and bb smaller and 
hence reduces the detection rate. Therefore, these new annihilation channels could 
either increase or decrease the detection rate, depending which of these two effects 
wins over. We found that, for small values of tan/3 and /i < 0, the WW channel 
can become dominant if open, basically because in this case the neutralino contains 
a rather large neutral wino component. This explains the distortion of the detection 
rate curves in the \x < half of Figs. ^ and |6]. The effect of the ZZ channel turns 
out to be negligible in the SU(5) x U(l) model, since neutralinos with m x > Mz 
have very small higgsino components. The same argument applies to the hA channel 
which sometimes opens up in the dilaton case for rather large values of tan/?. 

In the dilaton case, for large values of tan/3 the CP-odd Higgs boson A can be 
rather light, and the presence of the A-pole when m x ~ \rriA makes the relic density 
very small. £l x h\ as a function of m x , is first lower than 0.05, it increases with 
m x , and eventually reaches values above 0.05, when neutralinos move away from the 
A-pole. Thus, the capture and detection rates also show this behavior, which can 
be seen as the few "anomalous" lines in Figs. |5J and For lower values of tan/?, 
fi x /ig < 0.05, and there is no such effect. In the minimal SU(5) model, since the 
allowed points include different supersymmetry breaking scenarios and several values 
of mt and tan /3, all these features are blurred. Nonetheless, in the same range of m x 
and for same values of tan/3 and m t , we have found the results of these two models 
comparable, with the rates in the £77(5) x U(l) model slightly smaller due to the 
smaller relic density. 

It is clear that at present the experimental constraints from the "neutrino tele- 
scopes" on the parameter space of the two supergravity models are quite weak. In 
fact, only the Kamiokande upper bound from the Earth can be used to exclude re- 
gions of the parameter space with m x w m-pe for both models, in particular for the 
SU(5) x U(l) model, due to the enhancement effect discussed above. However, it 
is our belief that the results presented in this paper will be quite useful in the fu- 
ture, when improved sensitivity in underground muon detection rates become avail- 
able. An improvement in experimental sensitivity by a factor of two should be easily 
possible when MACRO |T7| goes into operation, while a ten-fold improvement is en- 
visaged when Super-Kamiokande IT3] announces its results sometime by the end of 
the decade. More dramatic improvements in the sensitivity (by a factor of 20-100) 
may be expected from DUMAND and AMANDA ||19|| , currently under construction. 
In addition, as recently argued [35], we think that perhaps the full parameter space 
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of a large class of supergravity models, including the two specific ones considered 
here, may only be convincingly probed by a detector with an effective area of 1 km 2 . 
It is interesting to note that, with a sensitivity improvement by a factor of 100, a 
large portion of the \x < half parameter space of the minimal ££7(5) model can 
be probed. Unfortunately, the remaining portion, with fluxes below ~ 10~ 4 , can 
hardly be explored by underground experiments in the foreseeable future. For the 
SU (5) x £7(1) supergravity models, in Figs. [7| and | we have plotted the allowed points 



in the (m x ±, tan (3) space. These points are those obtained originally in Refs. p9|, |3l|| , 
such that the neutrino telescope constraint is also satisfied; no other constraints have 
been imposed. It can be seen here that the small voids of points for m + ^100 GeV 
and a variety of values of tan/3 are excluded by the constraint from the "neutrino 
telescopes" . In these figures we have marked by crosses the points in the {m x ± , tan j3) 
plane that MACRO should be able to probe (assuming an increase in the sensitivity 
by a factor of two) for the no-scale and dilaton scenarios respectively. Finally, Figs. |] 
and [H] show that, using this indirect technique, Super-Kamiokande can cover nearly 
half of the parameter space of the SU(5) x £7(1) model, assuming that an improve- 
ment by a factor of about twelve can be achieved. As expected, the constraints from 
future "neutrino telescopes" will be strictest for large values of tan j3. 

To keep the above results in perspective, we now discuss some effects that have 
not been taken into account in our analysis. As we have shown, for the two super- 
gravity models, the capture of the neutralinos by the Earth is more important than 
by the Sun. In the Earth case, we have considered only the primary direct capture 
of the neutralinos from the galactic halo, in which a neutralino is trapped by the 
Earth's gravitational field only when its velocity falls below the escape velocity at a 
point inside the Earth, as a consequence of its interaction with the nuclei around. 
However, there exists yet another mechanism by which neutralinos can be captured 
by the Earth, namely, the secondary indirect capture, first studied by Gould 0. In 
this mechanism, a neutralino first loses only enough energy to be bound in a solar 
orbit ("orbit capture"), but not enough to be directly captured by the Earth. The 
orbit-captured neutralinos further weakly interact with the nuclei in the Earth, and 
a fraction of them subsequently can be indirectly captured into the Earth's core. It 
was found by Gould that, for Dirac neutrinos of mass 10-80 GeV, direct and indirect 
capture by the Earth are of the same order of magnitude, and the kinematic enhance- 
ment of the total capture rate becomes enlarged and broadened 0. We believe that 
this result should also apply to the case of neutralino capture, which means that the 
Earth capture rate that we have computed in this work may have been significantly 
underestimated. Although to our knowledge the indirect capture was also not consid- 
ered in recent works that dealt with Earth capture [jn], 13, T%\, we feel that this 
important mechanism should be taken into account in future analyses of this type. In 
fact, from Figs. |7| and || we see that even before MACRO, interesting constraints may 
already be extracted from the current Kamiokande upper limits, had we included the 
indirect capture in our calculation. 

Recently, after the calculations in this paper had been completed, we received 
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Refs. |37| in which some other issues of particle physics relevant to our work have 
been addressed. Given the fact that the neutralino in the two supergravity models we 
study is mostly a "mixed" state, in particular for m x ~ m Fe , according to Fig. 7 of 
Ref. []36|| , we expect that the results of Ref. |3H|| , when they become appreciable, could 
lead to a shift of our detection rates in both directions by about 10% or less. We have 
not redone our calculations of the capture rate using the cross section of Ref. |37j . The 
inclusion of the new neutralino-gluon scatterings [|37|, when they become appreciable, 
would increase the chance that neutralinos can be captured by both the Sun and the 
Earth, barring the normally small interference effects which sometimes could render 
the total cross section of Ref. [37] smaller than what it would be without these new 
scatterings f38j. Therefore, the effect of the results of Ref. |37j on the detection 
rate, could somehow enhance or compensate that of Ref. | 36fl , depending whether the 
latter is an increase or decrease of the detection rate. However, even as a conservative 
estimate, we do not expect that the combined effect of Refs. [37|] would alter our 
results quantitatively by more than 20%, a change not large enough to affect our 
conclusions qualitatively. 

We have also looked into the possibility that the high energy neutrino flux may 

m 



in the Sun. This issue has been addressed in a 



be altered by MSW oscillations 

recent paper jK]. From their results we conclude that for the neutrino energy range 



relevant to this work, to v e oscillations are negligible. The flux may also be 
altered by to v T vacuum oscillations during passage from the Sun to the Earth. 
This effect will be significant only if the v^-fr mixing angle is large. Large mixing is 
disfavored by the general GUT based see-saw arguments An analysis based on 
these considerations coupled with phenomenological arguments in the context of the 
SU(5) x U(l) model |^2j supports a value for sin 2 28 « 10~ 4 , which is much too small 
to affect the flux values predicted here. 



6 Conclusions 

We have explored the possibility of detecting supersymmetry indirectly through the 
measurement of upwardly-moving muons in underground detectors or "neutrino tele- 
scopes". These muons originate from high-energy neutrino interactions in the rock 
below the detector, and these neutrinos result from annihilation of neutralinos in the 
Earth or Sun cores. The latter would have been captured by these heavenly bodies 
if the neutralinos constitute an important part of the galactic halo — an important 
assumption which should not be overlooked. 

The present day experimental upper bounds on the muon flux are only weakly 
constraining. In fact, this is mostly because the large possible fluxes for light neu- 
tralinos (see for example Fig. [| for small m x ) have already been ruled out by the 
LEP lower bounds on the neutralino mass. Nonetheless, there is a region of parame- 
ter space with m x ± « 100 GeV (corresponding to m x w mFe) and a range of values of 
tan/?, which is excluded at the 90% C.L. However, expected increases in experimental 
sensitivity in the next few years should turn this technique into a very efficient way 
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of probing the parameter space of the specific supergravity models considered here, 
even surpassing the reach of traditional direct detection collider experiments. 

We should remark that even though our explicit computations apply only to the 
SU(5) and 577(5) x U(l) models, the correlations among sparticle and Higgs boson 
masses, which play such a fundamental role in the quantitative results, are common 
to a large class of supergravity models with radiative electroweak symmetry breaking 



15| . In this respect, the specific models considered here include values of £o = 
m o/ m i/2 which are small (£ = 0, no-scale), moderate (£ = 1/Vc$> dilaton), and 
large (£o 3> 1, minimal SU(5)), and therefore span the whole allowed range. Thus, 
our results: (i) represent a significant sampling of what would be obtained using an 
arbitrary selection of parameters in a generic supergravity model, and (ii) will allow 
to select the correct supergravity model when experimental data start to restrict the 
parameter space. 
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Figure Captions 



Figure 1: The neutralino capture rate for the Sun and Earth as a function of the 
neutralino mass in the minimal SU (5) supergravity model. 

Figure 2: The upwardly- moving muon flux in underground detectors originating from 
neutralino annihilation in the Sun and Earth, as a function of the neutralino mass 
in the minimal £77(5) supergravity model. The dashed lines represent the present 
Kamiokande 90% C.L. experimental upper limits. 

Figure 3: The neutralino capture rate for the Sun and Earth as a function of the 
neutralino mass in the no-scale ££7(5) x U(l) supergravity model. The representative 
value of m t = 150 GeV has been used. Note the depletion of neutralinos in the halo 
near the Z-resonance, and the enhancement in the Earth capture rate near the iron 
nucleus mass (52.0 GeV). 

Figure 4: The upwardly-moving muon flux in underground detectors originating 
from neutralino annihilation in the Sun and Earth, as a function of the neutralino 
mass in the no-scale SU(5) x U(l) supergravity model. The representative value of 
m t = 150 GeV has been used. The dashed lines represent the present Kamiokande 
90% C.L. experimental upper limits. 

Figure 5: Same as Fig. 3 but for the dilaton SU(5) x U(l) supergravity model. 

Figure 6: Same as Fig. 4 but for the dilaton SU(5) x U(l) supergravity model. 

Figure 7: The allowed parameter space of the no-scale SU(5) x £7(1) supergravity 
model (in the (m x ±,tan/3) plane) after the present "neutrino telescopes" (NT) con- 
straint has been applied. Two values of m t (130,150 GeV) have been chosen. The 
crosses denote those points which could be probed with an increase in sensitivity by 
a factor of two. 

Figure 8: Same as Fig. 7 but the crosses now denote points which could be probed 
with an increase in sensitivity by a factor of 12. 

Figure 9: Same as Fig. 7 but for the dilaton SU(5) x £7(1) supergravity model. 

Figure 10: Same as Fig. 8 but for the dilaton ££7(5) x £7(1) supergravity model. 
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